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Abstract 

This paper considers a multiple input multiple output (MIMO) two-way relay channel, where two 
nodes want to exchange data with each other using multiple relays. An iterative algorithm is proposed 
to achieve the optimal achievable rate region, when each relay employs an ampUfy and forward (AF) 
strategy. The iterative algorithm solves a power minimization problem at every step, subject to minimum 
signal-to-interference-and-noise ratio constraints, which is non-convex, however, for which the Karush 
Kuhn Tuker conditions are sufficient for optimaUty. The optimal AF strategy assumes global channel 
state information (CSI) at each relay. To simplify the CSI requirements, a simple amplify and forward 
strategy, called dual channel matching, is also proposed, that requires only local channel state information, 
and whose achievable rate region is close to that of the optimal AF strategy. In the asymptotic regime 
of large number of relays, we show that the achievable rate region of the dual channel matching and 
an upper bound differ by only a constant term and establish the capacity scaling law of the two-way 
relay channel. Relay strategies achieving optimal diversity-multiplexing tradeoff are also considered with 
a single relay node. A compress and forward strategy is shown to be optimal for achieving diversity 
multiplexing tradeoff for the full-duplex case, in general, and for the half-duplex case in some cases. 

TMs work was funded by DARPA through IT-MANET grant no. W911NF-07-1-0028. 
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Fig. 1. Two way relay channel communication protocol 

I. Introduction 

We consider a multiple antenna two-way relay channel as shown in Fig. [T| where two nodes Ti and 
T2 want to exchange information with each other with the help of a relay node and all the nodes are 
equipped with one or more than one antenna. The two-way relay channel models the communication 
scenario where the destination terminal also has some data to send to source terminal e.g. downlink 
and uplink in cellular communication, or packet acknowledgments in a wireless network. The general 
discrete memory less two-way relay channel was introduced in [1], and the multiple antenna two-way 
relay channel in [2]. In the literature, the two-way relay channel is also known by several other names, 
including the: bidirectional relay channel [3]-[5] and analog network coding [6]. 

A specific embodiment of a multiple antenna two-way relay channel that assumes half-duplex relays 
and the absence of a direct path between source and destination was proposed in [2]. An illustration is 
provided in Fig. [T] As shown in Fig. [T] in phase 1 or the first time slot, both terminals Ti and T2 are 
scheduled to transmit simultaneously while the relay receives. In phase 2 or the second time slot, the relay 
is scheduled to transmit while terminals Ti and T2 receive. The key idea with the two-way relay channel 
is that each terminal can cancel the interference (generated by its own transmission) from the signal it 
receives from the relay to recover the transmission from the other terminal. The idea is reminiscent of 
work in network coding [7], though note that here the coding is done in the analog domain, [6] rather 
than in digital domain [7]. In this paper we only consider multiple antenna two-way relay channel and 
for brevity, drop the prefix multiple antenna from here onwards. 

There has been a growing interest in finding the capacity region of the two-way relay channel with a 
single relay node [3]-[5], [8]-[13]. Achievable sum rate expressions (sum of the rates achievable from 
Ti T2 and T2 Ti links) have been derived in [2] and [3]-[5], for the half-duplex two-way relay 
channel, using amplify and forward (AF), decode and forward (DF) and compress and forward (CF) at 
the relay. It is shown that in a two-way relay channel, it is possible to remove the ^ rate loss factor in 
spectral efficiency due to the half duplex assumption on the nodes. For a general full-duplex two-way 
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relay channel with a single relay node (Ti, T2 and relay can transmit and receive at the same time) 
achievable rate regions are derived in [8] for AF, DF, and CP. For the AWGN two-way relay channel 
(no fading), using nested lattice coding and DF at the relay, the achievable rate region has been shown 
to be very close to the upper bound for all SNRs [10], [11]. Using the deterministic channel approach, 
the achievable rate region has been shown to be at most three bits away from the upper bound for the 
full-duplex two-way relay channel [9]. The capacity region of the two-way relay channel has also been 
studied in [12], [13], where in [13], it has been shown that in the low SNR regime the upper bound can 
be achieved by choosing a suitable relay mapping function, together with LDPC codes. The achievable 
rate region [3]-[5], [8]-[13] does not meet the upper bound [14], in general. Consequently, the problem 
of finding the capacity region of the two-way relay channel is currently open. 

The problem of finding the capacity region of the two-way relay channel becomes even more chal- 
lenging when there are multiple relay nodes that can help Ti and T2, and to the best of our knowledge 
has not been addressed in the literature. The problem becomes hard, because it is known that for the 
one-way relay channel with multiple relay nodes, DF does not work well [15], while the partial DF 
and distributed CP [15] lead to complicated achievable rate regions that are very hard to compute and 
analyze. The same conclusion holds true for the two-way relay channel; the only simple strategy that 
is well suited for multiple relay nodes is AF. With this motivation, in this paper we attempt to find the 
optimal relay beamformers that maximize the achievable rate region of the two-way relay channel with 
AF. For the one-way relay channel with multiple relays, optimal relay beamformers have been found 
[16], however, they are not known for the two-way relay channel. 

For the case when both Ti and T2 have a single antenna, and each relay has an arbitrary number of 
antennas, we solve the problem of finding optimal relay beamformers by recasting it as an iterative power 
minimization algorithm. The iterative algorithm, at each step, solves a power minimization problem with 
minimum signal-to-interference-noise (SINR) constraints, for which satisfying the Karush Kuhn Tucker 
(KKT) conditions [17], [18] are sufficient for optimality. We consider both the sum power constraint 
across relays, as well as an individual relay power constraint. The optimal AF solution requires each 
relay to have channel state information (CSI) for all relays and leads to an achievable rate region that 
cannot be expressed in closed form. 

For the case when each relay knows its own CSI, finding the optimal AF strategy is quite hard and 
intractable, even for the one-way relay channel case [16]. To remove the global CSI requirement, and 
to obtain a simple achievable rate region expression, next, we propose a simple AF strategy, called 
dual channel matching strategy, which works for any number of antennas at Ti and T2. In dual channel 
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matching, relay k transmits the received signal multiplied with (G^Hj^+H^*G^*), if the channel between 
Ti and relay k is Hj., between relay k and T2 is G^, between T2 and relay k is G^ and between relay 
k and Ti is H^. Using dual channel matching, we lower bound the achievable rate region of the optimal 
AF strategy, which is unknown for more than one antenna at Ti and T2, and bound the gap between 
the optimal AF strategy and the upper bound. The dual channel matching is quite simple to implement 
and its achievable rate region can be shown to be quite close (by simulation) to the optimal AF strategy, 
when Ti and T2 each have single antenna. 

We upper bound the capacity region of the two-way relay channel using the cut- set bound [19] on the 
broadcast cut Ti (T2), and ri, r2, . . . , vk, and the multiple access cut ri, r2, . . . , and T2 (Ti), over all 
possible two phase protocols (with different time allocation between first and second phase). We show 
that the gap between the upper and lower bound (dual channel matching) is quite small for small values 
of K. In the limit K ^ 00, v/e show that the gap is constant with increasing K, and thus estabhsh the 
scaling law [20] of the capacity region of the two-way relay channel, which shows that ^ log K bits can 
be transmitted from both Ti T2 and T2 — > Ti, simultaneously. 

We also consider the problem of finding relay transmission strategies to achieve the optimal diversity 
multiplexing (DM)-tradeoff [21] of the two-way relay channel with a single relay node, in the presence of 
a direct path between Ti and T2. The DM-tradeoff captures the maximum rate of fall of error probability 
with signal to noise ratio (SNR), when rate of transmission is increased as rlogSNR. The DM-tradeoff 
for the two-way relay channel is a two-dimensional region spanned by the (di2(^i2,^2i),c?2i('"i2,'"2i))> 
where di2 and ^21 are the negatives of the exponent of the probability of error from Ti T2 and 
T2 Ti, respectively, when Ti is transmitting at rate ri2logSNR and T2 at r2ilogSNR. The DM- 
tradeoff for the one-way relay channel has been studied in [22]-[25], where notably in [25], it has been 
shown that the CF strategy achieves the DM-tradeoff for both the full-duplex as well as the half-duplex 
case. The DM-tradeoff of the two-way relay channel has been recently studied in [26], where upper and 
lower bounds are obtained on the DM-tradeoff which are shown to match for the case when each node 
has a single antenna. 

We first consider the full-duplex two-way relay channel and show that a shghtly modified version of 
the CF strategy [27] achieves the optimal DM-tradeoff. More importantly, we show that di2{ri2,r2i) 
(c?2i('"i2,'"2i)) does not depend on r2i (ru) and the two-way relay channel can be decoupled into two 
one-way relay channels using the CF strategy. Then we consider the more interesting case of half-duplex 
nodes, where the achievable rate regions are protocol dependent. For the two-way relay channel it is not 
known which protocol achieves the highest possible rates [3]-[5]. We use a three phase protocol, where 
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in phase one Ti transmits to both the relay and T2, in phase two T2 transmits to both the relay and Ti 
and in phase three the relay transmits to Ti and T2. This three phase protocol makes use all the direct 
links between different nodes in a two-way relay channel. For this three phase protocol, we propose a 
modified CF strategy and show that it can achieve the optimal DM-tradeoff in some cases. We conjecture 
that our strategy can also achieve the optimal DM-tradeoff in general, but we are yet to prove it. 

Notation: The following notation is used in this paper. The superscripts -'^j* represent the transpose 
and transpose conjugate. M denotes a matrix, m a vector and mi the i*^ element of m. For a matrix 
M = [mi m2 . . . m„] by vec(M) we mean [m^ m^^ . . . m^]^. det{A) and tr(A) denotes the 
determinant and trace of matrix A, respectively. E denotes the expectation. || • || denotes the usual 
Euclidean norm of a vector and | • | denotes the absolute value of a scalar. is a m x m identity matrix. 
\X\ is the cardinality of set X. We use the usual notation for u{x) = 0{v{x)) if |^^| remains bounded, 
as X — > 00. X ~ CM{0,a) means x is a circularly symmetric complex Gaussian random variable with 
zero mean and variance a and x\y ~ CM{0,a) means given y, x is a circularly symmetric complex 
Gaussian random variable with zero mean and variance a. C^'^^ denotes the set of M x matrices with 

w.p 1 

complex entries. x„ > y denotes that the sequence of random variables x„ converge to a random 

variable y with probability 1. We use a = 6 to denote equality with probability 1 i.e. Prob.{a = 6) = 1 

w.p.l 

and < is defined similarly, /(x; y) denotes the mutual information between x and y and h{x) the 

w.p.l 

differential entropy of x [19]. To define a variable we use the symbol :=. 

Organization: The rest of the paper is organized as follows. In Section [nj we describe the two-way 



relay channel system model, the protocol under consideration and the key assumptions. In Section III 



we 



obtain the optimal AF strategy to maximize the achievable rate region of the two-way relay channel. In 
Section [IV| we introduce a simple AF strategy, dual channel matching, and lower bound the achievable 
rate region of the optimal AF strategy of Section III In Section [Vj we derive an upper bound on the 
capacity of the two-way relay channel capacity and compare it with the achievable rate region of the 



optimal AF strategy and dual channel matching. In Section VI we show that the CF strategy can achieve 
the optimal DM-tradeoff for full-duplex two-way relay channel, in general, and in some cases for the 



half-duplex case. Final conclusions are made in Section VII 



II. System and Channel Model 

In this section we describe the two-way relay channel system model under consideration, and then 
present the relevant signal and channel models. 
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Fig. 2. Two-way relay channel system model with two phase communication 



A. System Model 

For the first part of the paper Section III IV and [Vj we consider a wireless network where there 
are two terminals Ti and T2 who want to exchange information via K relays, as shown in Fig. |2] The 
K relays do not have any data of their own and only help Ti and T2 communicate. The K relays are 
assumed to be located randomly and independently so that the channel coefficients between each relay 
and Ti and T2 are independent. We also assume that there is no direct path between Ti and T2 and that 
they can communicate only through the K relays. This is a realistic assumption when relaying is used 
for coverage improvement in cellular systems, since at the cell edge the signal to noise ratio is extremely 
low for the direct path. In ad-hoc networks, it can be the case that two terminals want to communicate, 
but are out of each other's transmission range. 

We assume that both the terminals Ti and T2 have M antennas and all the K relays have N antennas 
each. We further assume that both the terminals and all the relays can operate only in half-duplex mode 
(cannot transmit and receive at the same time). The communication protocol is summarized as follows 
[2]. In any given time slot, for the first a fraction of time, called the transmit phase, both Ti and T2 are 
scheduled to transmit and all the relays receive a superposition of the signals transmitted from Ti and 
T2. In the rest (1 — a) fraction of the time slot, called the receive phase, all the relays are scheduled 
to transmit simultaneously and both the terminals receive. Both Ti and T2 are assumed to have power 
constraint of P, while for relays we assume two different power constraints, the sum power constraint 



7 



M 



y 




Fig. 3. Channel model for the two-way relay channel between Ti, T2 and relay k 




Fig. 4. System Model of the two-way relay channel with direct path for Section [Vl] 



where the sum of the power of all relays is < Pr or the individual power constraint where each relay 
has power constraint of Pr. 



For the second part of the paper, Section VI we assume a two-way relay channel with a single relay 
node and the presence of a direct path between Ti and T2 as shown in Fig. |4] We assume that Ti has 
mi antennas, T2 has m2 antennas, and the relay node has nir antennas. 



B. Channel and Signal Model 

Throughout this paper we assume that all the channels are frequency flat slow fading block fading 
channels, where in a block of time duration Tc (called the coherence time), the channel coefficients 
remain constant and change independently from block to block. We assume that Tc is more that the 
duration of time slot used by Ti and T2 to communicate with each other as described before. As shown 
in Fig. [3} let the forward channel between Ti and the k^^ relay be = [hi^ h.2k ■ ■ ■ i^Aik] and the 
backward channel between k*^^ relay and Ti be = [h^^ • • • ^kM^ ■ Similarly let the forward 
channel between k*^^ relay and T2 be = [g^i gfc2 • • • Ekivi] and the backward channel between T2 



8 



and the k relay be = [g^^ g2f, ■ ■ ■ EMk\- ^'^^ Section VI where the direct path between Ti and T2 
is considered, the channel between Ti and T2 is denoted by H12 and in the reverse direction by HJ2- 
We assume that Hfc, G C^^*^^,H^, G^ G C*^^^,Hi2 G C™^^™SH^2 G C^^^^^^ with independent 
and identically distributed (i.i.d.) CJ\f{0, 1) entries. 

For the first part of the paper Section [Tllj |IV] and |Vj we consider the following signal model. The 
X 1 received signal at the k^'^ relay is given by 

Hp Hp 

'''' " V M^''^^ V M^^^^ "'^ '-^^ 

if xi and X2 are the M x 1 signals transmitted from Ti and T2 to be decoded at T2 and Ti respectively, 
with E{xi*xi} = E{x2*X2} = M, P is the power transmitted by Ti and T2, respectively. The noise is 
the X 1 spatio-temporal white complex Gaussian noise independent across relays with E(nfcn^) = Itv- 
Relay k processes its incoming signal to transmit a iV x 1 signal = W^r/j with X^^j^ E{t^t/s} < Pr 
(sum power constraint) or Ejt^tfc} < Pr (individual power constraint) in the receive phase. The M x 1 
received signals yi and y2 at terminal Ti and T2, respectively, in the receive phase, are given by 

K 

yi = J]]H^tfc + zi, (2) 

k=l 
K 

y2 = ^Gfctfe + Z2, (3) 

k=l 

where zi and Z2 are M x 1 spatio-temporal white complex Gaussian noise vectors with E(ziz*) = 
E(Z2Z*) = Im. 

Throughout this paper we assume that both Ti and T2 perfectly know {H^, H^, Gfc, G^} \/ k, k = 
1,2, , ... K in the receive mode. To be precise, in the receive phase (i.e. when Ti and T2 receive signal 
from all the relays), Ti and T2 both know {Hfc,Gfc} and {H^,G^} V /c, = 1, 2, , . . . K. We also 
assume that no transmit CSI is available at Ti and T2, i.e. in the transmit phase Ti and T2 have no 
information about what the realization of and Gk is going to be when it transmits its signal to all 
the relays in the transmit phase, respectively. 

In this paper we assume different CSI assumptions at the relay. For finding the optimal AF strategy 



(Section IIIi we assume that each relay knows Hfc, G^, Gk, for all fe = 1, 2 . . . , K. To reduce the 



CSI requirements next, we present a simple AF strategy in Section IV where we assume that relay k 



only knows Hfc,G^,Gfc,H^. In Section VI we assume that the relay knows Hi, G^, Gi, H^', as well 
as H12, the channel coefficient between Ti and T2. 
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III. Optimal AF strategy for two-way relay channel 

In this section we will find optimal relay beamformers that maximize the achievable rate region of 
the two-way relay channel with AF, when Ti and T2 have a single antenna each, M = 1. For simplicity 
of exposition, in this section we consider the case when each relay nodes has a single antenna, = 1. 
Generalizations to > 1 are straightforward, and will be described later. 

To start with, because of single antenna restriction, the channel between Ti and relay k is denoted by 
/ifc and between relay k and T2 denoted by g^. For the reverse direction the channel coefficients are the 
same as in forward direction but with an added superscript r, e.g. channel coefficient between relay k 
and Ti is denoted by h^.. With AF strategy, each relay node transmits the received signal multiplied with 
Wk to both Ti and T2. Thus, if xi and X2 is the transmitted signal from Ti and T2, respectively, then 
the received signal at Ti, yi, and T2, 2/2 is 



where n^, \/ k = 1, . . . , K is CM{0, 1) noise added at relay k and zi, and Z2 are CAA(0, 1) added at Ti 
and T2. Since xi and X2 are known at Ti and T2, respectively, their contribution can be removed from 
the received signal at Ti and T2, respectively. Let the rate of transmission from Ti to T2 be R12 and 
from T2 to Ti be R21, then from Q 



Thus, the achievable rate region for the two-way relay channel with AF for a sum power constraint 
across all relays, i.e. pR = -P Ef=i(kfc^fcP + kfcS'I.P) + Ef=i kfcP < Rr is the set TZ^P.Pr) = 
^pr<Pr(.Ri2, R21) and for individual power constraint at each relay, i.e. p^R = P{\'Wkhk\'^ + \wkgl\'^) + 
|u^fcP < Pb. is the set Tl{P, Pr) = ^p^aKPa, k=i,...,K{Ri2, R2i)- Therefore, the problem is to find optimal 
that achieve the boundary points of the region Tl{P, Pr), for both the sum power constraint and an 
individual power constraint. 

For the one-way relay channel, no communication from T2 to Ti, optimal w^s have been found in [16] 
to maximize Ri2. The solution of [16], provides an upper bound on individual rates R12 and R21 and 



K 

VPhlwkglx2 + VPhlwkhkXi + hlwkUk + zi 

k=l 



K 

VPgkWkhkXi + \fPgkWkglxi + gkWknk + Z2 

k=l 



(4) 




V 
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is equivalent to solutions where R12 or R21 is greedily maximized disregarding the other. The problem 
in the two-way relay channel case is to find optimal w^'s such that Rgum = R12 + ^21 is maximized, 
for each /3 G [0, 1], where R12 = jSRsum, and R21 = (1 — f3)Rsum- Towards that end, we use the rate 
profile method [28] to identify Wk's that meet the boundary point of Tl{P,Pji). Next, we only consider 
the sum power constraint across the relays. For individual power constraints the same procedure can be 
applied as pointed out later. Thus, the optimization problem can be formulated as follows. 

Maximize^,,, k=i,2,...,K Rs 



I'Sum 



subject to log (^1 + U^K^j,,^,,.^ J > PRs 



PZkUl^khkl' + \wkgl\') + Ef=i \nJk\' < Pr. 
An equivalent problem to this problem is the following iterative power minimization problem subject 
to rate constraints, 

Minimize^^^ fc=i,2...,x PR = -PEf=i(kfc^fcP + kfcSfcP) + Ef=i kfcP 

subject to log [I + ' J > (6) 

where at each iteration R^um is changed to maximize the achievable rate, subject to power constraint. 
To be precise, if the value of -R"^^ at iteration i is say x and the solution to Q is feasible (i.e. if 
VR ^ Pr) then x is incremented in next iteration, otherwise decreased. Choice of the step size of 
increase or decrease determines the speed of convergence to the optimal rate R^^m^ for which < Pr. 
One possible starting point for R^um i^ ^ times the maximum R12 provided by [16] for one way relay 
channel. The step size can be chosen by bisection between the last feasible R^um (initially 0) and the last 
infeasible Rgum- Even though this equivalent problem provides a solution to ([5]) in a iterative manner, 
the problem (|6]) is in general non-convex, and not easy to solve. To overcome this limitation, we recast 
the problem Q as a standard power minimization problem subject to signal-to-interference-noise ratio 
(SINR) [18], where the forwarded noise from each relay plays the role of interference. For a given f3 



'For an individual power constraint the same can be done by checking at each iteration whether the obtained solution pr is 
feasible with individual power constraints or not. 
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and Rgum' the problem ([6]) is of the form 
Minimize 

subject to 



k=i,2-,K PR = PIlk=ii\wkhk\ + kfcSfcl ) + Efc=i kfcl 



> 
> 



P 

2(l-/3)H^„„ 



70 

-1 



(7) 



This problem again is non-convex, however, it is of the form 



TI- 



CS) 



Minimize f{x) 

subject to ||ai(x)|p - < 0, V i, 

where f{x) is a convex function, ai{x) is an affine function of x and bi{x) > V i, by noting that if 
SfcLi 9kWkhk, or J2k=i ^k'^kdl are less than zero or complex, then they can be scaled by appropriate 
phases to make them real and positive, without changing the objective function or the constraints |^ 

For the problem ([8]l, it has been shown in [18], that if the problem is strictly feasible, then KKT 
conditions [17] are necessary and sufficient to find the optimal solution. It is easy to see that the problem 
(|7]) is strictly feasible and therefore KKT conditions are sufficient for optimality. The Lagrangian of 
problem ([7]) is of the form 

1 



where w = [wi 



wAw* + Ai ^wBw* 
. . wk] and 



cw 



T|2 



70 



+ 1 + A2 wUw 



lew^P + 1 



+ 1 











|2 I |„r|2\ 



B 



and c 



51^1 



\9i\ 





gxhK 








P{\h2V + \glV) + l 




71 







\9k\ 



D 














... 























Differentiating the Lagrangian yields 

A + AiB + A2D 



— c c H e e | w 

70 71 



0, 



^An immediate consequence of this property is that the optimal solution does not change if all g'f.s are scaled by e-''^^, or all 
hl.'s are scaled by e^*^^. 
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and the optimal w is found by solving for Ai and A2 using the constraints |j 

Therefore, by recasting our original problem of obtaining the boundary points of Tl{P, Pr) to the 
power minimization problem with SINR constraints, we have shown that the optimal solution can be 
found in an efficient way. In Section |V} we plot the achievable rate region of the optimal AF strategy 
and compare it with the lower bound obtained by using dual channel matching, and an upper bound. 

Recall that we only considered a two-way relay channel, where each relay had a single antenna, = 1. 
Extension to > 1, is straightforward by replacing gkWkhk by gfcWfchfe, gkWk by gfcWfc, hlwkQk by 
h^Wfcgfc and /i^Wfc by hJ^W^, which are scalars as before, and the optimal solution to W^'s can be 
found using the iterative power minimization algorithm ([6]). 

Our algorithm to optimize the achievable region with AF is fairly simple, however, it assumes that 
each relay has CSI for all the relay nodes, and requires M = 1. Finding optimal relay beamformers 
where each relay has only its CSI, and M > 1, is rather complicated and has not been solved even for 
the one-way relay channel [16]. Another limitation of the optimal AF strategy is that the expression for 
the obtained rate region cannot be written down in close form, and therefore does not allow analytical 
tractability for comparison with an upper bound. To remove these restrictions, in the next section we 
propose a simple AF strategy, called dual channel matching, where each relay uses its own CSI, and for 
which the achievable rate region expression can be written down in a closed form. Since dual channel 
matching is in general, a suboptimal AF strategy, the achievable rate region of dual channel matching 
lower bounds the rate region of the optimal AF strategy, and allows to estimate the difference between 
the optimal AF strategy and the upper bound. 

IV. Dual Channel Matching Strategy 

In this section we propose a simple AF strategy, called dual channel matching, and derive a lower 
bound on the achievable rate region for the two-way relay channel. With the dual channel matching 
strategy relay k multipUes ^//3fc(G^H^ + H^*G5^*) to the received signal and forwards it to Ti and T2, 
where fik is the normalization constant to satisfy the power constraint. Dual channel matching tries to 
match both the channels which the data streams from Ti to T2 and T2 to Ti experience at each relay 
node. The motivation for this strategy is that for one-way relay channel (i.e. T2 has no data for Ti) with 
one relay node, the optimal AF strategy is to multiply V2DU^ to the signal at the relay, where the 
singular value decomposition of Hi is UiDiV* and Gi is U2D2V2 and D is a diagonal matrix whose 
entries are chosen by waterfilling [29]. In dual channel matching the complex conjugates of the channels 

^Clearly, the optimal w lies in the null space of some matrix that is a function of A, B, c, e, and D and hence not unique. 
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are used directly rather than the unitary matrices from the SVD of the channels [29]. This modification 
makes it easier to analyze the achievable rates for the two-way relay channel. Note that the dual channel 
matching is an extension of the listen and transmit strategy of [30] for the one-way relay channel, where 
each relay transmits the received signal after scaling it with the complex conjugates of the forward and 
backward channel coefficients. 

Together with dual channel matching we restrict the signal transmitted from Ti and T2, xi and X2, 
respectively, to be circularly symmetric complex Gaussian distributed with covariance matrix Ejxixi*} = 
E{x2X2*} = Im> to obtain a lower bound on the achievable rate region of two-way relay channel. 
Moreover, we use a = | i.e. Ti and T2 transmit and receive for same amount of time. The achievable 
rates R12 and R21 using the dual channel matching can be computed as follows. 

From ([1]), the received signal at the k^^ relay is given by 

Hp rp 
" V M^^^^ ^ V M^^^^ ^ 

Using dual channel matching as described above, at relay k, G^H^ + H^*G^* is multiplied to the 
received signal so that the transmitted signal is given by 

tfc = -xZ/Jfe (G^Hfc + H^*G") rfc 

where (5^ is to ensure that Y^^=i ^*k^k = Pr |^ With dual channel matching the received signal at T2 is 
given by 

K 



y2 



^Gfctfc+z. (9) 



k=l 



Expanding Q we can write 



k=l k=l 

^ V ' 

A 

K 

+ V%Gk {GlUl + Ul*Gl*) n, + z. 

A— 1 

Since X2 and all the channel coefficients are known at T2, the second term can be removed from the 
received signal at T2. Moreover, as described before xi is circularly symmetric complex Gaussian vector 
with covariance matrix Q = Im, thus the achievable rate for Ti to T2 link is [31] 



R12 = ^/(xi; y2) = \ log det (im + AA* B^B^ + 1^^ j , (10) 



''This is for the sum power constraint. For an individual power constraint, /3 is chosen such that t^t^ = P^j for each k. 
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since Ejnfcn^} = E{zz*} = Im, V k. Similarly, we obtain the expression for R21, 

i?2i = ^logdet I lAf + CC* I > 'DfcD^ + lAf 1 I, (11) 



^ log det (im + CC* DfcD^ + Im^ 



where C = Zt=i VtHI (G^H* + H^Gr) G^ and = (G*H* + H^GD • This rate 

region expression obtained is analytically tractable and can be used to compare the loss between the 
optimal AF strategy and the upper bound. Another interesting question of interest is how does the 
achievable rate region behaves with K. To answer that question, we turn to asymptotics and compute 
the rate region in the limit ii' ^ cx), in the next lemma. 
Lemma 1: As K grows large, K —>■ 00, 

M 

lim Ru =^ — logK + 0(l), 



M 

lim R21 = — logK + 0(l). 
w.p.l 2 



Proof: Consider 



2i?i2 - log det KIm = log det + AA* B^B^ + Ia/^ j - log det KIm, 



/ K 

1 A A* ^ 



-IN 



= log det I^-Ia/ + (^E^'^^fc 

To satisfy the sum power constraint, let /5 = |^ where ci is a constant such that 

ci = E {((G^H*, + HrGD + HrGD rfc} , 

which is same for all k. Then, 

K 



/ , {^k"-k + "-k ^k ) ^fc, 



^ V ciM K ^ 
which by using strong law of large numbers, converges to. 



A w.p.l I PPr 2t 

Im, 



VK V ciM 

since E{GfcG^} = E{H^Hfc} = NIm V k, and E{GfcH^*} = OIm V k. Same result holds true for 
— With B = 

^ Po 1 ^ 

BfeBfc - — — 2^ (Gfe (GfcH;!, + Hfc Gfc jj (Gfc (G^H^ + G^ jj , 

k=l ^ k=l 



^Equal power allocation among relays. 
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which again using the strong law of large numbers converges to ^6Im, for some finite 9, since 
Hfc, Gjt, HJI, G^B^ are i.i.d. with finite variance. Thus, in the limit K ^ oo, 

2Ri2 - logdet/^lM ^ Mlog ( X^'^^y , ) , 

\M{Pr9 + ci)J 

and thus it follows that 

Ri2 =^^logK + 0{l). (12) 

Similarly we get the achievable rate R21 on the T2 to Ti link as 

M 

lim R21 = —logK + 0{l). (13) 

A— >oo w.p.l Z 



Discussion: In this section we introduced the dual channel matching AF strategy, and obtained a lower 
bound on the capacity region of the two-way relay channel. Dual channel matching is a simple AF 
strategy that requires local CSI, and as we will see in Section [Vj has achievable rate region very close to 
that of the optimal AF strategy (Section [In]) for M = 1. We also derived the asymptotic achievable rate 
region of the dual channel matching, by taking the limit K 00, and using the law of large numbers. 
We showed, that in the asymptotic regime, both R12 and R21 scale as ^ log K with increasing K. 

Next, we derive an upper bound on the capacity region of the two-way relay channel, and compare it 
with the achievable rate region of the dual channel matching. 

V. Upper Bound on the Two- Way Relay Channel Capacity 

In this section we upper bound the capacity region of the two-way relay channel using the cut-set 
bound [19] for the broadcast cut, and the multiple access cut. We assume a general two-phase protocol 
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where for a fraction of the time slot Ti and T2 transmit to all relays and the rest of the (1 — a) fraction 
of time slot all relays simultaneously transmit to both Ti and T2. Note that to lower bound the capacity 
of the two-way relay channel using dual channel matching, we used a = ^ which might be suboptimal. 
We prove later that for the asymptotic case of K 00, a = ^ is optimal. 

The upper bound is derived as follows. We start by first separating Ti and then T2 from the network 
and apply the cut set bound [19] to upper bound the rate of information transfer between Ti T2 and 
T2 —>■ Ti, respectively. Using the cutset bound, we first show that the maximum rate at of information 
transfer from Ti T2 (T2 ^ Ti) is upper bounded by the maximum rate of information transfer between 
Ti (T2) and ri,r2, ■ ■ ■ ,rK (broadcast cut) and also by the maximum rate of information transfer between 
ri, r2, . . . , rx and T2 (Ti) (multiple access cut). Fig. [5] and [6] Then we use the capacity results from [31] 
to upper bound the maximum rate through the broadcast cut for the case when CSI is only available at the 
receiver (all relays) and all the relays collaborate to decode the information. Similarly, for the multiple 
access cut as shown in Fig. [6j we upper bound the maximum rate at which all the ri,r2, ■ ■ .r^ can 
communicate to T2 (Ti) by using capacity results from [31], when CSI is known both at the transmitter 
(all relays) and the receiver (Ti,T2) and all the relays collaborate to transmit the information. 

Broadcast cut - To derive an upper bound we make use of the cutset bound (Section 14.10 [19]). 
Separating the terminal Ti from the rest of the network and applying the cutset bound on the broadcast 
cut as shown in Fig. |5] 

R12 < a{/(xi;ri,r2,...,ri^,y2|ti,t2,...ti^,X2)}. (14) 

Again applying the cutset bound while separating the terminal T2, 

R21 < a{/(x2;ri,r2, . . . , r^, yi|ti, t2, . . .ti^,xi)} (15) 

for some joint distribution p(xi, ti, t2, . . . , tx, X2), where R12 and R21 are the maximum rates at which 
Ti can communicate to T2 and T2 can communicate to Ti respectively, reliably. By the definition of 
mutual information [19] 

/(xi;ri,r2, . . • , r^, y2|ti, t2, . . .ti<-,X2) = /(xi;ri,r2, . . . ,ri^|ti,t2, . . .ti<-,X2) 

+ /(xi;y2|ri,r2,...,rii-, ti, t2, . . . tj^, X2). (16) 

By expanding the mutual information in terms of entropy, 

/(xi;ri,r2,...,rx|ti,t2,...t/^,X2) = /i(xi|ti, t2, . . . ti^, X2) 

- /i(xi|ri,r2, . . . ,ri^,ti,t2, . ..tK,X2) 
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Since conditioning can only reduce entropy [19], 

/(xi;ri,r2,...,rii-|ti,t2,...tj^,X2) < /i(xi|x2) 

- /i(xi|ri,r2, . . . ,rA',ti,t2, . . .tx,X2). 
Note that ti, t2, . . . , ti^ is a function of ri, r2, . . . , rx, which implies 
/(xi;ri,r2, . . . ,rx|ti,t2, . . .t/^,X2) < /i(xi|x2) 

- /i(xi|ri,r2, . . . ,r/^,X2) 

and henc43 

I(xi;ri,r2, . . . ,ri^|ti,t2, . . .ti^,X2) < /(xi; ri, r2, . . . ,ri^|x2). 
Given perfect channel knowledge at terminal T2, 

/(xi;y2|ri,r2, . . • ,rx, ti,t2, . . .ti^,X2) = /(xi,Z2) 
where Z2 is the AWGN noise. Since xi and Z2 are independent, /(xi,Z2) = 0, and therefore from (16 

I(xi;ri,r2, . . . , r^^, y2|ti, t2, . . .tx,X2) < /(xi;ri,r2, . . . ,rA-|x2). 



(17) 



Hence from (14), 



R12 < /(xi;ri,r2, . . . ,rx|x2). 
Similarly, by interchanging the roles of xi and X2, 

R21 < /(x2;ri,r2, . . . ,rK|xi). 



(18) 



(19) 



Therefore it is clear that both R12 and i?2i is upper bounded by the maximum information flow through 
the broadcast cut Fig. [5] when all the relays are allowed to collaborate. Expanding the mutual information 
in terms of differential entropy, 

/(xi;ri,r2,...,rx|x2) = /i(ri, r2, . . . , r;^|x2) - /i(ri, r2, . . . , rj^ |xi, X2). 



From ([T]), 

Since is known at relay k, 



P P 



P P P \ 

(ri, r2, . . . , rx|x2) = ^ I V V M^^^^ + n2, • • • , y J^^^^^ + t^k\^2 1 • 



^Without X2, in [20], this inequality has been shown to be an equality, which is incorrect. 
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Since conditioning can only decrease entropy, 



/i(ri,r2, . . . ,rx|x2) < h 



Hixi + ni 



-H2X1 + n2, 



M V M 

With perfect knowledge of and at relay k, 

/i(ri,r2, . . . ,ri^|xi,X2) = /i(ni,n2, . . ■ ,nK) , 

and it follows that 




HxXi + UK 



/(xi;ri,r2, . . . ,rx|x2) < h 



-H1X1 + ni 

M 

/i(ni,n2, . . .,nK) ■ 



P / P 

H2X1 + no, ■ ■ ■ , \ H;<-Xi + IlK 

M ' ' V M 



(20) 



Thus, we have shown that R12 is upper bounded by the maximum rate from Ti to ri, . . . , r/< without any 
interference from T2 and when all r^'s can collaborate to decode the message, which is quite intuitive. 
Using results from [31] when CSI is known only at the receiver, the R.H.S. of (20l is upper bounded by 
logdet (Im + Ef=i ffH^Hfe^, which implies 



/(xi;ri,r2, . . • ,rx|x2) < alogdet + ^ -^H^H^^ 



and therefore, from (fTSll 



R12 < a/(xi;ri,r2, . . . ,ri^|x2) < alogdet ( Im + ^ — H^H^ | 



Interchanging the roles of xi and X2 and replacing with G^, 



K 



R21 < Q;/(x2;ri,r2, . . . ,rx|xi) < alogdet Im + ^ 



fc=i 



KM 



k 



(21) 



(22) 



(23) 



Multiple access cut - Again by using the cutset bound, we bound the maximum rate of information 
transfer R12 {R21) from Ti —>■ T2 (Ti — > T2) by the maximum rate of information transfer across the 
multiple access cut as shown in Fig. [6| Using cutset bound, R12 and i?2i are bounded by 



Now, 



R12 < (1 - a)/(xi,ti,t2, . . . ,tx;y2|x2) 
R21 < (1 - a)/(x2,ti,t2, . . . ,tA-;yi|xi). 

/(Xi,ti,t2, . . . ,tx;y2|x2) = /l(y2|x2) - /l(y2|tl,t2, . . . ,t/^,X2) 

+ /l(y2|tl,t2, . . . ,tif,X2) - /l(y2|tl,t2, . . . ,tK,Xi,X2) 



(24) 
(25) 
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K Relay Terminals 



Fig. 6. Multiple Access Cut 



Note that given ti, t2, . . . , tx, y2 is independent of xi and X2, 

/i(y2|tl,t2, . . . ,t/^,Xi,X2) = /l(y2|tl,t2, . . . ,tx,X2) = /l(y2|tl,t2, . . • ,t;^). 

Therefore 

/(Xi,ti,t2, . . . ,tx;y2|x2) = /i(y2|x2) - /l(y2|tl,t2, . . . ,ti^). 
Since conditioning can only reduce entropy, 

I(xi,ti,t2, . . . ,tA';y2|x2) < h{y2) - /l(y2|tl,t2, . . . ,ti^), 
and by definition of mutual information 

I(xi,ti,t2, . . . ,ti^;y2|x2) < /(ti,t2, . . . ,ti^,y2). 



Hence from ( |24| ), 

Ri2 < (l-a)/(ti,t2,...,tK;y2). (26) 

Following similar steps we can also bound R21 as, 

R21 < (l-a)/(ti,t2,...,tx;yi). (27) 

Thus, R12, R21 are bounded by the maximum rate of information from ri, . . . , r^^ to Ti or T2. Next, we 
compute the maximum rate of information from ri, . . . , to Ti or T2. Recall from Q that the received 
signal y2 is given by 

K 

y2 = ^Gfctfc +Z2. 

k=l 

Note that 

y2 



I(ti, t2, . . . , ti^; y2) = / ( ti, t2, . . . , ti^-; 
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Dividing y2 by yK, we get 

This can also be written as 

y2 1 



y2 1 Z2 



[Gi G2 ... G;,][tit2...tK]^ + ^ 
is. y Sx V -IV 



Note that <I> is a M x NK matrix. Now assuming that all the relays know Gfc, VA; (allowing cooperation 
among all relays), with sum power available across all relays bounded by Pr, we have from [31], 



min {NKM} 

i?i2 <(!-«)/( ti,t2,...,ti^;^) < (1- a) Y max{0,log(i^Azi/)} (28) 



where A;, / = 1, 2, . . . , min {NK^ M} are the eigen values of <I><I>* matrix and v is chosen such that 

max{0,z^ - — } = Pr. 

— A/ 

i=\ ' 

Similarly, one can obtain the bound for i?2i by replacing G^ by H^. 



Combining ( p2| ), ( |23| ) and ( [28] ), gives the upper bound on the capacity region of the two-way relay 
channel. Comparing the upper bound with the lower bound obtained using the dual channel matching 
( 10|11| ), one can see that they do not match for any arbitrary value of K. In the asymptotic regime, 
however, they can be shown to be only an 0(1) term away as K ^ oo, as proved in the next Theorem. 
This asymptotic result implies two things, one that the performance of the dual channel matching, and 
consequently the optimal AF strategy (which we don't know for M > 1), does not degrade in comparison 
to the upper bound with increasing K, and two, it provides us with the capacity scaling law of the two-way 
relay channel. 

In Figs. |7] and [8j we plot the achievable rate region of the optimal AF strategy, the lower bound 
obtained using dual channel matching, and the upper bound for K = 2 and K = 4, with M = 1, = 1 
and P = Pr = lOdB with sum rate constraint across relays. Note that the achievable rate region of the 
optimal AF region is symmetric, as expected, because of the symmetry in parameters of communication 
in both directions in a two-way relay channel. Another important point to note here is that, the achievable 
rates of dual channel matching are quite close to that of the optimal AF strategy, even though it uses 
only local CSI. Thus, dual channel matching is a good candidate for AF in practical implementation of 
two-way relay channels. Also notice that the difference between the upper and lower bound is less than 
the 3 bit bound of [9]. 

Next, we prove that the lower bound (dual channel matching) and the upper bound on the achievable 
rate region of the two-way relay channel are only an 0(1) as K oo. We prove the theorem by 
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Comparison of optimal AF, Dual Channel Matching, and Upper Bound for M=1 , N=1 , K =2 



- Upper Bound 



-H — Optimal AF 
-O — Dual Channel Matching 




2 bits/sec/Hz 



Fig. 7. Comparison of upper and lower bound of the capacity region of the two-way relay channel with K = 2,M = 1^N = 1, 
P = Pn = IQdB 



Comparison of optimal AF, Dual Channel Matching, and Upper Bound for M=1 , N=1 , K =4 



Upper Bound 

H — Optimal AF 

© — Dual Channel Matching 




Fig. 8. Comparison of upper and lower bound of the capacity region of the two-way relay channel with K = A, M = 1, N = 1, 
P = Pji = lOdB 
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approximating the upper bound in the K oo and comparing it with the asymptotic lower bound 



obtained in ( 12 13 1. 



Theorem 1: The upper and lower bounds on the capacity region of the two-way relay channel differ 
by a 0(1) term as K ^ oo, and the capacity scaling law is given by 

Rl2 < yl0gK + 0(l), 
Rl2 < yl0gi^ + 0(l). 



Proof: We first approximate the broadcast cut upper bound (|22]| as ^ cx). From (22) 



Ri2 < a/(xi;ri,r2, . . . ,ri^|x2) < alogdet I Im + ^ ^H^H^ | . (29) 



Consider 



/ K p \ ( I I ^ P 

logdet [Im + Y. M^k^k - logdeti^lM = logdet -1m + J^Yl M^^^'^ 

V k=l / \ k=l 

Using strong law of large numbers 

k=l 



and it follows that 



logdet(ilM + i|:^H;H,]^Ml„g(^) 



which using ([22]) implies 



lim Ri2 < aM log K + 0{1), (30) 

X— >oo w.p.l 



since M, A^, P are finite integers. Similarly, 



lim R21 < aM log K + 0{1). (31) 

K^oo w.p.l 



Next, we approximate the upper bound of the multiple access cut. From (28), 

, , min{NK,M} 

i?i2 < (l-a)/(ti,t2,...,t;^;^) < (1-a) V max{0, log (A^^i^)}, (32) 
V VKJ 

where A;, / = 1, 2, . . . , min {NK, M} are the eigen values of matrix and u is chosen such that 

min {NK,M} ^ 

max{0, - — } = Pr. 
1=1 ' 



By definition = Ylk=i ^fc^fc- Frorn strong law of large numbers 

K 

K'-^oo K 



k=l 
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Therefore 

Xi = N Vi = l,2,...M, ^ ^ = + 1 



and from (|32]l 



and consequently, as — > oo 

Ri2 < (1 -a)M log K + 0{1), (33) 



w.p.l 



and similarly 



Combining ( |30|31| ) and ( |33p4 



i^ai < (1 -a)M log 7^ + 0(1). (34) 

ui.p.l 



M 

i?i2 < min{a,l -a}M log K + 0(1) < —log 7^ + 0(1), 



2 

M 



R21 < min{a,l -a}M log K + 0(1) < —log i^ + 0(l). (35) 



2 

Comparing ( [35] ) to the asymptotic lower bound (12 13 1 we conclude that (a) upper and lower bounds 
on the capacity region of the two-way relay channel differ by a 0(1) term as ^ 00, and (b) the 
capacity scaling law is given by 

M 

R12 < ylogK + 0(l), 
R21 < y logic + 0(1). 



To illustrate the result of Theorem [T] in Figj9} we compare the lower (dual channel matching) and 
upper bound on the sum rate R12 + R21, and show that they both scale similarly with increasing K for 
M = 2, = 1, P = Pji = lOdB with sum rate constraint across relays. 

Discussion: In this section we obtained upper bounds on the capacity region of the two-way relay 
channel, and compared it with the dual channel matching lower bound. To compute the upper bound we 
used the cut-set bound and the capacity results of [31]. The lower and upper bound expressions do not 
match in general, however, in the asymptotic case, where the number of relays are large, K ^ 00, we 
showed that they are only an 0(1) term away from each other. Thus, the dual channel matching and 
consequently, the optimal AF strategy are almost optimal in the asymptotic regime. For the finite number 
of relay nodes (finite K), we use Monte Carlo simulations to quantify the gap between the lower and 
the upper bound. From Figs. |7] and [8} we can see that gap between the lower (dual channel matching) 
and upper bound is rather small, and inside the 3 bit bound of [9]. 
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Sum rate comparison with increasing number of relay nodes K, M=2, N=1 

14 
12 
10 

5 8 

o 
w 

~w 

S 6 
4 
2 


5 10 15 20 25 

Nr of relays K 

Fig. 9. Capacity scaling of two-way relay channel with AI — 2, N = 1, P = Pr = lOdB. 

Another important observation to make is that the lower bound with dual channel matching was obtained 
using a = 5 i.e. Ti and T2 transmit and receive for equal amount of time. Since this lower bound is only 
a 0(1) term away from the upper bound (arbitrary a), distributing equal amount of time for transmit 
and receive phase is optimal in achieving the right capacity scaling. 

Compared to the asymptotic results on the one-way relay channel [20], [32], our results show that 
by two-way relay channel one can remove the ^ rate loss factor on the capacity, which comes from 
the half-duplex assumption on the terminals and relays. Therefore with two-way relay channel one can 
achieve unidirectional full-duplex performance with half-duplex terminals. 

VI. Diversity-Multiplexing Tradeoff 

In this section we consider a two-way relay channel with a single relay node, and characterize its 
DM-tradeoff. We consider both the full-duplex and half-duplex nodes, where Ti and T2 have nii and 
m2 antennas, respectively, and the relay node has rrir antennas. An important difference in this section 
from the previous ones is the presence of direct link between Ti and T2 as shown in Fig. |4] 

To characterize the DM-tradeoff, for both the full-duplex and half-duplex case, we first obtain an upper 
bound on the DM-tradeoff and then propose a modified CF strategy to achieve the upper bound. We first 
discuss the full-duplex case followed by the half-duplex case. 




H Lower Bound Dual Channel Matching 

■© — Upper Bound 
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A. DM-tradeoff of Full-Duplex Two-Way Relay Channel 

The signal model for this section is as follows. Let xi, X2 and x,. be the signal transmitted from Ti, 
T2 and the relay node, respectively. Similarly, Let yi, y2 and y^. be the signal received at Ti, T2 and 
the relay node, respectively. Recall that channel coefficient between Ti and relay node is denoted by H, 
between Ti and T2 is denoted by H12, between the relay node and T2 is denoted by G, where note 
that, compared to previous sections, we have dropped the subscript index of relay node, since we only 
consider one relay. All the channel coefficients in the reverse direction (right to left) are denoted by 
channel coefficient in the forward direction (left to right) with an added superscript r, e.g. the channel 
coefficients between the relay node and Ti is denoted by H''. Let the transmit power at Ti, T2 and the 
relay node be Then, 

yi = \ — H12X2 + W — H''xr+ni, 

V m2 V nT'r 

y2 = \ — H12X1 + W — Gxr + n2, 

V mi V TTT-r 

yr = \ — Hxi + W — G''x2 + nr. (36) 

V mi V ^2 

Let the rate of transmission from Ti to T2 and T2 to Ti be R12 and R21, respectively. Following [21], 
let Ci2(SNR) and C2i(SNR) be the family of codes, one for each SNR for transmission from Ti to 
T2, and T2 to Ti, respectively. Then we define ri2 (r2i similarly) as the multiplexing gain of Ci2(SNR) 
if the data rate i?i2(SNR) ( i?2i(SNR)) of Ci2(SNR) (C2i(SNR)) scales as ri2 (r2i) with respect to 
log SNR, i.e. 

i?i2(SNR) 

lim ,^ = ri2 

SNR-»oo log SNR 

and di2{ri2,r2i) ((i2i(?'i2, ?'2i)) as the rate of fall of probability of error Pei2 (^e2i) of Ci2(SNR) 
(C2i(SNR)) with respect to SNR, i.e. 

Pei2(SNR) =SNR-'^^^(''^^ 

The negative of the SNR exponent of the error probability ^12(^12, r2i) or d2i{ri2,r2i) captures the DM- 
tradeoff, where di2{ri2,r2i) {d2i{ri2,r2i)) is the maximum diversity gain possible from Ti to T2 (T2 to 
Ti) for a given ri2 and r2i. Note that the error probability Pei2(SNR) and Pe2i(SNR) are functions of 
both ri2 and r2i because of simultaneous transmission between Ti and T2. 

Next, we upper bound the DM-tradeoff of the two-way relay channel, the region spanned by (^12(^12, ^21) 
and d2i{ri2,r2i), by allowing cooperation between Ti and relay, and T2 and relay node. 



^Having different transmit power constraints for Ti, T2 and the relay node do not cfiange the DM-tradeoff. 
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Lemma 2: The DM-tradeoff of a two-way relay channel is upper bounded by 

di2{ri2,r2i) < min{(mi - ri2)(m,. + 7712 - ri2), (mi + - ri2)(m2 - ri2)}, 

d2i{ri2,r2i) < min{ (7712 - r2i) (777,. + 7771 - r2i), ("1-2 + 777^ - r2i)(r77i - r2i)}, Vri2, r2i. 
Proof: We will prove the lemma only for (ii2(?'i2, ?'2i)- For d2i(7'i2, ?^2i) it follows similarly. Consider 
the case when T2 has no data to send to Ti. This assumption can only improve di2(ri2, r2i). Then 
first assume that the relay node and T2 are co-located and can cooperate perfectly. In this case, the 
communication model from Ti to T2 is a point to point MIMO channel with 7771 transmit antennas 
and nir + m,2 receive antennas. The DM-tradeoff of this MIMO channel is (7771 — ri2){mr + 7772 — 
ri2), and since this point to point MIMO channel is better than our original two-way relay channel, 
di2{ri2,r2i) < {mi — ri2){mr + 7n,2 — ri2) . Next, we assume that Ti is co-located with relay node 
and both of them can perfectly cooperate for transmission to T2. This setting is equivalent to a MIMO 
channel with mi + m,. transmit and 7772 receive antenna with DM-tradeoff (7771 + 777,. — ri2)(7772 — ri2). 
Again, this point to point MIMO channel is better than our original two-way relay channel and hence 
di2('^i2, ?'2i) < ("7i + 777^ — ?'i2)(m,2 — ri2), which Completes the proof. ■ 
To achieve this upper bound we consider the CF strategy [27], with a slight modification and prove 
that it is sufficient, to achieve the optimal DM-tradeoff. We make few changes to the original CF strategy 
[27] to suit the two-way relay channel communication, which are as follows. Let the rate of transmission 
from Ti to T2 and T2 to Ti be R12 and R21, respectively. Instead of generating only one codebook 
at Ti as in [27], both Ti and T2 generate 2"^i" and 2"^"^ independent and identically distributed 
and X2 according to distribution p{x^) = YYi=iPi^^'d ^^'^ Pi^"^) = nr=i^'(^2i)' respectively. The 
codebook generation at the relay and the relay compression and transmission remains the same as in [27], 
i.e. the relay node generates 2"^° independent and identically distributed x^ according to distribution 
l?(x") = nr=i?'(^") label them Xr{s), s G [1,2"^°], and for each Xr{s) generates 2""^ y's, each 
with probability p{y\xr{s)) = Y[7=iP(y''-\^ri{s)). Label these y{z\s),s G [1,2"^"] and z G [1,2"^] and 
randomly partition the set [1,2"^] into 2^^° cells Sg, s G [1,2"^°]. Let in block i the message to send 
from Ti is Wi, and from T2 is Vi, then Ti sends xi{wi), T2 sends X2{vi) and the relay sends Xr{si) if 
Zi G Si, where y{zi\si-i),yr{i — 1), Xr{si-i) are jointly typical. Decoding at both Ti and T2 remains 
the same as in [27], however, note that in this case Ti knows xi{wi) and T2 knows X2{vi) apriori and 

^This upper bound is valid as long as the coherence time Tc is smaller than the time it takes for T2 to compute the channel 
coefficients and feed them back to Ti, which is at least mi + m2 [33]. Otherwise, T2 can help Ti in acquiring transmit CSI, 
for which case, potentially infinite diversity gain can be achieved [34], violating the present upper bound. 
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therefore can use them to decode Vi and Wi respectively. This strategy has been previously considered in 
[8] to obtain achievable rate region. 
With this two-way CF strategy, the following rates are achievable, 



Ri2 < /(xi;y2y|xrX2), 
R21 < /(x2;yiy|xrXi), 

with the compression rate constraint 

max{/(yr; ylx^xiyi), /(y^; y|xrX2y2)} < min{/(xr; yi|xi), /(x^; y2|xrX2)}. 



(37) 



(38) 



The rate region and the compression constraint are a little different from [27]. The rate region differs 
because of conditioning by xi or X2, which is due to the prior knowledge of xi at Ti and X2 at T2. The 
new compression rate constraint incorporates the condition that the quantized version of y,., y^ can be 



decoded at both Ti and T2. In the next Theorem we compute the outage exponents for (37i and show 
that they match with the exponents of the upper bound. 

Theorem 2: CF strategy achieves the DM-tradeoff upper bound (Lemma |2]|. 



Proof: To prove the Theorem we will compute the achievable DM-tradeoff of the CF strategy ( 37 1 and 
show that it matches with the upper bound. 



To compute the achievable rates subject to the compression rate constraints for the signal model (36 1, 
we fix y = yr + n^, where ng is rrir x 1 vector with covariance matrix Nlm^- Also, we choose xi, X2, 



and independent 



and Xr to be complex Gaussian with covariance matrices ^Imi, ^Im2> and 
of each other, respectively. Next, we compute the various mutual information expressions to derive the 
achievable DM-tradeoff of the CF strategy. By the definition of the mutual information 



/(xi;y2y|xrX2) 



/l(y2y|XrX2) - /l(y2y|XrX2Xi). 



From ( [36| ), arranging y2 y in a vectorized form we get 

^-Hi2Xi + 



y2 




y 





-GXr + n2, 



Hxi + J ^G'"x2 + + ng. 



and consequently 



/i(y2y 1x^X2) = log-L'i^, 



(39) 



(40) 



where 



= det I — * + 

mi 



N + l^ I„„ 
I 



m2 



and = [H12 H]. 
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Moreover, from f39\ 



/i(y2y|xrX2Xi) = logdet 



(iV + 1)1^^ 
I™. 



which implies 



Similarly, one can show. 



I(xi;y2y|xrX2) = log 



I(x2;yiy|xrXi) = log 



where 



L^^ = det 



P 

mi 



I, 



(iV + 1)™- 

r rl 

(iV + 1)™- 




(41) 



and H^2' = [H^2 Gl- 



(42) 



Next, we compute the value of that satisfies the compression rate constraints (38 1. By the definition 
of mutual information. 



/(yr.;y|xr.X2y2) = /i(y|xrX2y2) - /i(y|x,.X2y2yr) 

= /i(yy2|XrX2) - /l(y2|XrX2) - /l(y |XrX2y2yr) ■ 



(43) 



From (40l, /i(yy2|xrX2) = logL^^^. From signal model (36l, it is easy to see that /i(y2|xrX2) = 
logLi2, where L12 = det ^^Hi2H*2 + Ima) ■ Given y^, y has only the noise term n^, and hence 
/i(y|xrXiyiyr) = logiV™'-. Therefore, from (43), 



/(yr;y|xrX2y2) = log 



Similarly one can compute 



I(yr.;y|x.rXiyi) = log - — §- — , where L21 = det — H12H12 + Irm 



m2 



(44) 



(45) 



Again using the definition of mutual information, 

I(xr;yi|xi) = /i(yi|xi) - /i(yi|xrXi) 

= log^L - logL21, 



(46) 



where = det + ^H^H" + j, since yi = ^4H^2X2 + ^£H^x, + m. Simi 

larly. 



/(xr;y2|x2) = log 



r2 

Ll2 



(47) 



where = det (^Hi2H*2 + £GG* + I, 
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To satisfy the compression rate constraints (38 1, from (44l, ( |45] ), (46 1, ( [47] ), clearly 

LI' \ 



N > 



mill I log g^, log g^} 



(48) 



We choose N to satisfy the equality (48 1. From [21], to compute (ii2(n2> ?'2i)> it is sufficient to find 
the negative of the exponent of the SNR of outage probability at T2, where outage probability at T2, 
^oMi(''i2logSNR), is defined as 



From (37 41 1, 



R 



12 



log 



P{Ri2 < ri2logSNR) 



where N is given in (48 1. Then, 



Po^.t(ri2logSNR) 



P lo 



P 



(7V+ l)m 



(49) 



Choose / G Z such that {N + 1)"^- < I 



(iV + 1) 

l/rUr 



< SNR*^^^ 



< ri2 log SNR 



in (48 1. Then, 



Po«t(ri2logSNR) < P 



+ 1 ) , where N is such that it meets the equaUty 



+ 1 



< SNR' 



P 



P 



(Lf)V-^(L2jiA 



/I/™. ((Lf )Vrn. + (Lfjl/'n.) 



< SNR'^^ 



< P 



(^r2)l/m.(^2jl/m„ 
(^r2)l/m,-(^2jl/m„ 



< ;i/™,-SNR''i^/™'- 



< 



(50) 

(51) 
(52) 
(53) 



(Lf)iM- + (Lf,)i/ 

where the last equality follows because multiplying SNR by a constant does not change DM-tradeoff. 
From here on we follow [25] to compute the exponent of the Po^f (ri2 log SNR). Let 



/ P 



r2TTr2* 



Then clearly from (|40|, L^^ < L^^, therefore using Lemma 2 [25], it follows that 



-r2 



PoutirulogSm) < P 



(L-f)V--(Lf,)V-- 
(L^2)i/-. + (L2^)i/" 



< SNR*^^^/*^ 



(54) 



(55) 
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Moreover, notice that for non-negative random variables X and Y and a constant c [25], P[XY / [X + 
Y) <c)< P{X < 2c) + P{Y < 2c), thus, 

Pout{ri2 log SNR) < P ((L^?)^/™-- < 2SNR"^^/™-) + P ((L?,)^/™'- < 2SNR"^^/"^'-) , (56) 

= P {Lfi < SNR"^^^) + P {Ll^ < SNR'^^^) , (57) 

= SNR-'^^^''^^) + SNR"^^(''^^\ 

= SNR"™™''-'^'*-^'^-*' '^^(''12)1. (58) 

Therefore, to lower bound the DM-tradeoff we need to find out the outage exponents di{ri2) and ^2(^^12) 
of L^l and Lf^,. Notice that, however, L^f is the mutual information between Ti and T2 by choosing 
the covariance matrix to be ^Imi|^ and allowing the relay and T2 to cooperate perfectly. From [21], 
choice of i^^mi as the covariance matrix does not change the optimal DM-tradeoff, therefore, di{ri2) = 
{mi —ri2){mr + m2 — ^12). Similar argument holds for L\^., by noting that is the mutual information 
between Ti and T2 if the relay and Ti were co-located and could cooperate perfectly, while using 
covariance matrix ^^^^^ I^i+m^- Thus, d2{ri2) = {mi + nir - ri2){m2 - ru). Thus, for Ti to T2 
communication, the achievable DM-tradeoff with CF strategy meets the upper bound (Lemma [2]). A 
similar result can be obtained for T2 to Ti communication by choosing an appropriate n G Z such 



that {N + !)*"'■ < n ( ( ^ ) + 1 ) , where N is such that it meets the equality in (|48| and by 



carrying out the outage exponent analysis of R21 = log ^^^^^^^ and lower bounding L"^^ by L2], where 

B. Half-Duplex Two-Way Relay Channel 

In this section we compute the DM-tradeoff of the half-duplex two-way relay channel where all the 
nodes (Ti, T2 and the relay) are half- duplex. For the half-duplex case, the achievable rate regions are 
protocol dependent and the optimal protocol is unknown in general [3]-[5]. Here we compute the DM- 
tradeoff of a three phase protocol, that is intuitively optimal (difficult to prove), where for ti fraction of 
the time slot Ti transmits to both T2 and the relay, t2 fraction of the time slot T2 transmits to Ti and 
the relay, and for the rest (1 — ti + 12) fraction of the time slot the relay transmits to both Ti and T2. 

For this communication protocol the rates R12 and R21 are upper bounded by the following expressions. 

R12 < max mill {ii/(xi;yr,y2),ti/(xi;y2) + {I - ti - t2)/(xr;y2)} , 
R21 < maxmin{t2i'(x2;yr,yi),t2/(x2;yi) + (1 - ti - t2)/(xr;yi)} , 

tl ,t2 



taking the role of SNR. 
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where the first argument in the minimum is obtained by allowing the relay and the T2 (Ti) to collaborate in 
the receive mode, and the second argument is obtained by simply adding the maximum mutual information 
possible at T2 (Ti) while in receiving mode. Using the rate region expression, we define the upper bound 
on the DM-tradeoff of the half-duplex two-way relay channel as follows. 



From the definition of (54i, 



P(ti/(xi;y„y2) <ri2logSNR) = P (ti log L^^ < ri2 log SNR) , 

:= SNR-'^"('^^^\and (59) 
P (ti/(xi; ys) + (1 - ti - t2)/(x,; ys)) = P {h log L12 + {1 - h - ts) log Lr2 < n2 log SNR) , 

:= SNR-'^"-("^=\ (60) 
where L2r = det (l^^ + ^GG* j . Thus, di2(ri2,r2i) < maxt^^t, min {4c (^2), d^^acii" 12)} ■ Simi- 



larly we can obtain upper bound for ^21(^12,^21) by replacing ti by t2 in d59 60 1 



To achieve this upper bound we consider the CF strategy of subsection VI-A except that in this case 
the compression signal y is chosen such that it is jointly typical with the received signals y^ti and y^j^ 
received in time ti and t2 from Ti and T2, respectively With this CF strategy the achievable rate 
region is given by 

-R12 < ti/(xi;y2y|xr,X2)), 
-R21 < t2i'(x2;yiy|xr,xi)), 
subject to the following compression rate constraint 

(ti + t2)max{/(yr;y|xrXiyi),/(yr;y|xrX2y2)} < (1 - (^i + t2)) min{/(xr; yi|xi), /(x,.; y2|xrX2)}. (61) 

To compute these rates, we let xi, X2 and x,. to be the same as in the full-duplex case and y = 
Yrti + Yrta + whcrc iig is the complex Gaussian vector with zero mean and covariance matrix iVI,.. 



Following the same steps as in (48 1 to (57i, we obtain 



P{Ri2 < ri2 log SNR) < P{ti log Llf < ri2 log SNR) + 

^f {2{h+t2)-l)tl ^ , (l-(tl+t2))tl ^ ^2 ^ 1 qmr\ 

P — — log L12 H — — log Lf,. < ri2 log SNR , 

V *1 + 12 11+12 / 

:= SNR-'^"^^^^) + SNR-'^™-("^^\ 

= SNR"™''^'^^'^'''')''^""'^^'''')}. (62) 

'"in [4] a similar strategy has been proposed, but there, two separate compression signals are chosen that are jointly typical 
with yrti and yrt2 individually, and then a deterministic function of the two compression signals is transmitted from the relay, 
which results in a different rate region expression from the one obtained here. 
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Thus the achievable (ii2(ri2, ^-21) < maxf^^t^ min {ci^^(ri2), dm^cC'^ia)}- Note that the expression for 
'^12('^12, ?'2i) is independent of r2i, and because of symmetry in R12 and R21 expressions, similar bounds 
can be obtained for R21 by replacing ti with t2, and is given by 

P(i?2i < rai log SNR) < P(ti log L^j < r2i log SNR) + 

/(2(ti+i2)-l)t2i . ^ (l-(ti+i2))t2i .1 . 1 ciMpA 

P — — log L12 + — — log L^^ < r2i log SNR , 

\ ti +12 ti +12 / 

:= SNR-'^"(''=^^) + SNR-"^— ("^^^ 



which implies 

rf2i(n2,^^2i) < maxmm\dll{r2i),d'^'^^{r2i) \ . (63) 

It is clear that the lower bound (62 63 1 and the upper bound ( 59|60 ) on the DMT of the half-duplex 
two-way relay channel do not match for the general case. By comparing the achievable DM-tradeoff and 
the upper bound, the next Theorem characterizes the cases for which CF strategy is optimal. 

Theorem 3: The proposed CF strategy achieves the optimal DM-tradeoff of the half-duplex two way 
relay channel if 

• the bottleneck of the channel is the broadcast cut, i.e. d^^(ri2) < c^mac(^i2) and correspondingly in 
the upper bound (i^^(ri2) < c^mac(^i2)> and with similar relation for (i^^(r2i) and (i^jj^(r2i)also. 

. otherwise if (^^^ij^]^ log L12 + log Lf^ = tilogLi2 + {I - h - t2)Lr2, and with 

similar relation for T2 to Ti communication. 
Proof: Follows immediately by comparing the lower bound (62i and the upper bound ( 59|60 l on the 
DM-tradeoff. ■ 

Discussion: In this section we showed that the CF strategy achieves the optimal DM-tradeoff of the 
two-way relay channel for the full-duplex case, in general, and for the half-duplex case in some cases. For 
both the full-duplex and half-duplex case we upper bounded the DM-tradeoff allowing different nodes 
to collaborate with each other while transmitting or receiving. For the full-duplex case, we modified the 
CF strategy of [27] and showed that it decouples the two-way relay channel into two one-way relay 
channel and achieves optimal DM-tradeoff on each of the two one-way relay channels. For the half- 
duplex case, as observed before, the achievable rate region and consequently the DM-tradeoff depends 
on the communication protocol. We used a three phase protocol that makes use of all the direct links 
between Ti, T2, and the relay. For the three phase protocol we proposed a modified CF strategy where 
the compression signal is chosen such that it is jointly typical with the signals received at the relay node 



"The same strategy can also be found in [4] 
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in phase 1 and 2. Using this CF strategy, we obtained a lower bound on the DM-tradeoff that is shown 
to match with the upper bound under some conditions. For the general case also, we believe that the 
proposed CF should be optimal in terms of achieving the DM-tradeoff, however, showing that is quite 
difficult because of the different mutual information quantities involved as well as the maximization over 
the time durations of phase 1 and 2. 

Our result for the full-duplex case is similar to [25], where it is shown that the CF strategy achieves 
the optimal DM-tradeoff in one-way relay channel. For the half-duplex case, however, because of three 
phase communication protocol and added compression rate constraints we are unable to reach the same 
conclusion of [25] in general, that CF achieves the optimal DM-tradeoff in half-duplex one-way relay 
channel. 

VII. Conclusion 

In the first part of the paper, we addressed the problem of finding optimal relay beamformers to 
maximize the achievable rate region of the two-way relay channel with multiple relays, when each relay 
uses AF. The use of AF strategy is motivated by the fact that all the other known relay strategies such as 
DF, partial DF and CF, do not work well in the presence of multiple relays, and moreover, AF is quite 
simple to implement. 

For the case when both the terminals Ti and T2 have a single antenna and each relay has an arbitrary 
number of antennas, we found an iterative algorithm to compute the optimal relay beamformers. The 
algorithm is equivalent to solving a power minimization problem subject to SINR constraints at each 
step. The power minimization problem at each step is non-convex, however, for which it is sufficient to 
satisfy the KKT conditions to obtain the optimal solution. 

The derived optimal AF strategy maximizes the rate region with AF, but is restricted to the case of 
a single antenna at Ti and T2, and cannot be extended easily for the multi-antenna case. Moreover, it 
also requires each relay to have global CSI, and does not have a closed form achievable rate region 
expression. To relax the single antenna restriction and global CSI requirement, we then proposed a dual 
channel matching strategy, which requires local CSI, and showed that the gap between the rate region of 
the optimal AF and dual channel matching is quite small when both Ti and T2 have a single antenna. The 
dual channel matching works for any number of antennas at Ti and T2, and has a closed form expression 
for the achievable rate region. We then compared the achievable rate region of the dual channel matching 
with an upper bound to quantify the loss while using dual channel matching. The analytical expressions 
of the lower and the upper bound did not match, and we used simulations to show that the gap is quite 
small. In the asymptotic regime of — >^ 00, however, using the analytical expressions, we proved that 
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the achievable rate region of the dual channel matching, is only a constant term away from the upper 
bound. Thus, we obtained the capacity scaling law for the two-way relay channel. Compared to [20], 
[32], our capacity scaling law for the two-way relay channel shows that with two-way relay channel, 
there is a two-fold increase in the capacity compared to unidirectional communication. 

In the second part of the paper, we considered the problem of finding coding strategies that achieve the 
optimal DM-tradeoff in a two-way relay channel with a single relay node, in the presence of direct path 
between Ti and T2. We showed that the CF strategy achieves the optimal DM-tradeoff of the full-duplex 
two-way relay channel, by first decoupling the two-way relay channel into two one-way relay channels, 
and achieving the optimal DM-tradeoff on each of the two one-way relay channel. For the half-duplex 
case we showed that a modified CF strategy for a three phase transmission protocol achieves the optimal 
DM-tradeoff for some cases. 
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